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ABSTRACT 

We compare the efficiency of weak lensing-selected galaxy clusters counts and of the weak lensing 
bispectrum at capturing non-Gaussian features in the dark matter distribution. We use the halo model 
to compute the weak lensing power spectrum, the bispectrum and the expected number of detected 
clusters, and derive constraints on cosmological parameters for a large, low systematic weak lensing 
survey, by focusing on the fim-crg plane and on the dark energy equation of state. We separate the 
power spectrum into the resolved and the unresolved parts of the data, the resolved part being defined 
as detected clusters, and the unresolved part as the rest of the field. We consider four kinds of clusters 
counts, taking into account different amount of information : signal-to-noise ratio peak counts; counts 
as a function of clusters' mass; counts as a function of clusters' redshift; and counts as a function of 
clusters' mass and redshift. We show that when combined with the power spectrum, those four kinds 
of counts provide similar constraints, thus allowing one to perform the most direct counts, signal- 
to-noise peaks counts, and get percent level constraints on cosmological parameters. We show that 
the weak lensing bispectrum gives constraints comparable to those given by the power spectrum and 
captures non-Gaussian features as well as clusters counts, its combination with the power spectrum 
giving errors on cosmological parameters that are similar to, if not marginally smaller than, those 
obtained when combining the power spectrum with cluster counts. We finally note that in order to 
reach its potential, the weak lensing bispectrum must be computed using all triangle configurations, 
as equilateral triangles alone do not provide useful information. The appendices summarize the halo 
model, and the way the power spectrum and bispectrum are computed in this framework. 
Subject headings: cosmological parameters - large-scale structures - weak gravitational lensing 



1. INTRODUCTION 

Since its first detec ti ons i n the early 2000's 
(iBacon. Refregier fc EUisI l2000t Ivan Waerbeke et all 
120001 : IWittman et al.ll200d|) . weak gravitational lensing, 
the coherent distortion of distant galaxies by inter- 
vening dark matter, has become a premier tool to 
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and has been sho wn to be the most p r omising probe 
of d ark energy (|Albrecht et al.l 120061 : iPeacock et al.l 
12006*). Surveys' optimization and systematics mini- 
mizations in both software and hardw are haye been 
investigated (iHevmans et al . 2006; Mass ev et al.l l2007al: 
Amara & Refregierll2007', '2008; Paulin-Henri ksson et al.l 
2008; Amara ct al. 2009; Bridle et al. 2009), favoring 
wide surveys, with well-controlled, stable Point Spread 
Function, with comprehensive photometric redshift 
follow-up. Those characteristics are shared by ambitious 
upcoming large area surveys, such as the Large Synoptic 
Survey Telescope (LSST)^, the Panoramic Survey 
Telescope & Rapid Response System (Pan-STARRS)*^, 
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Euclid'' and the Joint Dark Energy Mission (JDEM)®. 
They will provide us with a large amount of high quality 
imaging, well fitted to measure dark energy with weak 
lensing. Anticipating a significant improvement in re- 
moving systematics, we are left with the question of how 
to best extract the cosmological information out of the 
data. For instance, one must decide how to optimally 
capture non-Gaussianities and break the degeneracies 
between cosmological parameters as constrained by 
the extensively studied and measured weak lensing 
power spectrum. In this paper, we ignore primordial 
non-Gaussianities, and consider non-Gaussianities due 
to the growth of structures only. 

The power spectrum, the Fourier transform of the 2- 
point correlation function of the shear field, has been 
the most used measuremen t so far, both from g round 
(e.g. iMassev et al .! 12005"; IVan Wae rbeke e t all 120051: 
Jarvis et al.ll2006f" Hoekstra et al.l 2006: Scmb oloni et al.l 
20061; iBenjamin et al. 2007; Fu et al. 2008) and from 
space (e.g. Rhodes ct al. 2004; Hc ymans et al. 20051; 
ISchrabl3ack et al.ir2007l : iMassev et al.ii20'07b( ) . The intro- 
duction of tomography, the three-dimensional informa- 
tion of the shear field, that captures structure evolution, 
has tightened constraints on cosmological parameters, in 
particular the matter density Om and the amplitude of 
density fluctuations erg (|Massev et al.ll2007bl ). Despite its 
success, the power spectrum leaves us with well known 
degeneracies between parameters, that we must break by 

http:/ /sci.esa.int/science-e/ www/area/index.cfm?fareaid=102. 
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combining it with other measurements and/or probes. 
The power spectrum capturing only the Gaussian fea- 
tures of the field, it is natural to introduce measures 
of non-Gaussianities, and combine them with the power 
spectrum, to break degeneracies and tighten constraints. 

Higher-order correlations measure the matter den- 
sity field's non-Gaussian features. The lowest one, 
the weak lensing 3-point correlation function (3PCF), 
and its Fourier transform, the bispectrum, have been 
given a lot of attention in th e past few years 
(ISchncidc r fc Lombardi|[200 3: Takad a fc Jahi 2003a,b^d; 
Zaldarriaga fc Scoccimarro' 200 3: .Takada fc Jain. i2QqJ; 
Schneider ct al. 2005; Bcnab ed fc: Scoccirnarrol 120061: 
Semboloni et al.ll2008t i Joachimi et al.ll2009l : IVafaei et al.l 
20091) . However, a clean measurement of the 3PCF is 
difficult, thus few paper s have reported observa tional 
measuremen t s so far. iBernardeau et al.l (|2003( ) and 
iJarvis et al.l ()2004( ) have measured the skewness, but a 
full measurement of the shear 3PCFs, or of the conver- 
gence bispectrum, is still to be done. Several efforts are 
underway, including an algorithm to measure the conver- 
gence bispectrum directly in Fourier space (jPires et al.l 
|2009a). 

Large-scale structures such as clusters of galaxies, 
the result of the non-linear evolution of density fluc- 
tuations, are the non-Gaussian features that we want 
to take into account. Assessing their abundance as 
a function of various parameters, such as redshift, 
has been known as a powerful probe and used as 
this for several years (e.g. Oukbir fc Blanchard 199^ 
Eke et all '1998; W ang fc Steinhardti [l998t iHolder et al 
2001 1 : jRefregier et al.l 12001 iBattye & Welle: 
2003[ IPierpaoli et al.l l2003t IW ang ct al. 2004; 
Horellou fc Be rgd 120051: iMarian fc Bernstein i200 



Gladdcrs ct al. 120071: iMantz et al.r T2008: Sahlc n et alj 
20091). In particular, tight constraints can be obtained 
flm and (78 ■ Because the weak lensing clusters 
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selection function is rapidly evolving with redshift, 
counting weak-lensing-selected clusters is less sensitive 
than using catalogs of X-ray, optical, or Sunyaev- 
Zel'dovich (SZ) clusters. Despite this fact, constraining 
cosmology is possible with weak- l ensing clusters only 
(IWeinberg fc Kamionkowskl l2003t iMarian fc BernsteinI 
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studied by sev e ral authors Ce.g. |Cooray fc Hul l2001al : 
iTakada fc JahU [200l iTakada fc Bridlel 120071 ). Most 
constraints' predictions have so far been made using 
a fitting func tion for the non-li n ear power spectrum, 
such as Peacock fc Doddsl (| 19961 ): iMa fc Frvl ^2000) or 
ISmith et al.l ()2003l ). However, fitting functions tuned 
for a ACDM uni verse must be used w ith caution when 
one IS varying w (|Joudaki et al.ll2009D . We use the halo 
model to compute the weak lensing statistics. Despite 
its simplicity, it is well suited to cosmological parame- 
ters prospects, since it allows us to vary the dark energy 
without extrapolating fitting functions out oi w = — 1. 
Appendix lAl summarizes the halo model description that 
we use in this paper. The halo model allows us to de- 
fine the power spectrum as the separated contributions 
of the resolved and unresolved parts of the weak lensing 
field. Such a separation, since it gets rid of highly non- 
Gaussian features, can help improve constraints on cos- 
mological parameters. We consider four different types 
of cluster counts, with more or less intrinsic cosmolog- 
ical information : counts as a function of shear signal- 
to-noise ratio (S/N) only, counts as a function of mass, 
counts as a function of redshift, and counts as a func- 
tion of mass and redshit. We show that when com- 
bined with the power spectrum, counting clusters just 
as a function of their S/N gives constraints similar to 
those obtained when combining the power spectrum with 
counts of clusters taking the full mass and redshift in- 
formation. We then show how the weak lensing bis- 
pectrum captures non-Gaussian features and breaks the 
power spectrum degeneracies as well as clusters counts. 
We assume a Euclid-like survey, 20,000 deg^ wide, with 
rig = 40 galaxies arcmin"^ with median redshift = 1, 
without any external (e.g. from CMB) priors, so as to 
show what can be done with a weak lensing survey alone. 

Section [2] summarizes the methods we employ to con- 
strain cosmological parameters with different weak lens- 
ing statistics. Section |3| presents our results. We con- 
clude in section|4l Our halo model c ode will be made part 
of the public icosmo package^ (Refr egier et al.ll2008| ). 

2. METHOD 

2.1. Weak lensing selection function 

The weak lensing selection function for clusters of 



Marian et al.| [2009; Mat uri et al. | |2009|; jWang et al.| galaxies has already been extensively studied, by using 



20091) ■ as shown on real data by iBerge et al.l ()2008l ) 



Although such counts give weaker constraints than 
the weak lensing power spectrum, combining them 
with the p ower spectrum efficiently breaks degeneracies 
(jPires et a l. 2009b). 

In this paper, we show how combining the weak lens- 
ing power spectrum with the weak lensing bispectrum or 
counts of weak-lensing-selected clusters of galaxies will 
allow one to tighten constraints on cosmological param- 
eters with upcoming weak lensing surveys. In particular, 
we show how we will be able to measure the equation of 
state w = P/ p, where P is the pressure, and p the den- 
sity of the dark energy, down to the percent level. While 
other authors studied the combi nation of t he power spec- 
trum with cluster counts Ce.g . ITakada fc Bridle 2007) 
and with the bispectrum (e.g. [Takada & Jain 2004), we 
present, for the first time, a detailed and consistent com- 
parison of both combinations at once. 

Constraints on cosmological parameters have been 
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selection function, based on a matched-filter approach 
for the signal-to-noise ratio created by a halo. We ne- 
glect projection effects and intrinsic alig nments , which 
have been investiga t ed e.g . bv Maturi et al. I ([2005), 
IMarian fc BernsteinI (|2006( ) and|Fan| (|2007t ). That selec- 
tion function has been introduced, without the mathe- 
m atical justification t hat follows, and tested on real data 
in IBerge et all (|2008[ ). 

We consider a spherically symmetric weak lens, decou- 
pled from its surrounding and alone along the line of 
sight, characterized by its convergence Kobs(^), where 
is the distance from its center (the convergence can be 
replaced by the shear without any loss of generality). 

^ http://www.icosmo.org 
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Fig. 1. — Weak lensing selection function for clusters of galaxies, 
in the redshift - mass plane. We assume a survey with Ug = 40 
galaxies per arcmin'^ , distributed along equation |4] Contours de- 
note the S/N of galaxy clusters. 



The number density of background sources is Ug. We as- 
sume that the noise originating from Poisson and 
intrinsic shape noise only, averages to in circular shells 
and that its variance in the circular shell i is given by 
< nf >= oi^jNi, where is the r.m.s error on the shape 
measurement and iV, is the number of lensed galaxies in 
shell i. We fit a theoretical model k(G) to the observable 

^obs 

The signal-to-noise ratio created by the lens is defined 
as the ratio of an estimator of to its associated er- 
ror, V = {Kv^ jcj (^(^Kw^^- The estimator is defined by 

(^Kw^ — J d^6'Kobs(^)w(0), where w{9) is a weight func- 
tion to be adjusted so as to optimize the S/N. This is the 
case when w = n, giving the optimized S/N as 




(1) 



which is consistent with previously published expres- 
sions. 

We will consider NFW halos (jNavarro et al.ll996( ) only, 
the S /N of which is expressed as (appendix [B]) 



j,y V^.^. (2) 

where G(c) = dx xg{x)^ is a function of the halo's 
concentration c only, well fitted by 

G(c) « ^ - ^ + 0.388 - 5 X lO-^c - 2.8 x 10" V 
c c 

(3) 

for 1 ^ c ^ 200, and the function g{x) is defined by 
Eq. dHlJ. The quantity (Z) in Eq. @ describes the 
effect of the distribution of galaxy sources. In the above 
equations, is the angular-diameter distance to the 
lens, Ps and parametrize the NFW halo (see Appendix 
[X)) . and Scrit,oo is defined in Appendix IB] 



Fig. 2. — Weak lensing convergence power spectrum, when one 
single redshift bin is used to measure it. The black curve shows 
the power spectrum when the contribution of all halos is taken 
into account. Other curves make use of halos whose signal-to-noise 
is less than a given threshold (the unresolved part) : i^tjj = 3 
(dashed), i-'tj, = 7 (dash-dot) and = 11 (dash-dot-dot-dot). 



We para metrize the redsh ift distribution of source 
galaxies by (jSmail et al.lll994 



n{z) = z°' exp 



(4) 



with a = 2, /? = 1.5, and zq ~ Zm/1.412, where Zm is the 
median redshift of the survey. 

Figure [T] shows our fiducial survey's selection function, 
in the mass-redshift plane. Contours originate from the 
S/N of halos in that plane. 

2.2. Weak lensing statistics 

2.2.1. Power spectrum 

The halo model allows us to decompose the full three- 
dimensional matter power spectrum along 

Ps{k) = pr{k)+pr[k), (5) 

where P™^''{k) is computed on the unresolved part of the 
halos distribution, defined as the ensemble of halos with 
S/N smaller than a certain threshold vtii\ Ps'^^ is com- 
puted on the resolved part of the halos distribution, the 
ensemble of halos with v > vth ■ Those names come nat- 
urally from the separation made on a mass map between 
peaks (defined as resolved structures with v > i^th) and 
the rest, unresolved part of the map. The resolved part 
can be seen as the most non-Gaussian features of the 
field. 

The full weak lensing power spectrum Pk.(£) follows 
the same decomposition. Figure shows the full power 
spectrum together with power spectra computed on un- 
resolved parts defined by the thresholds fth = 3, fth = 7 
and vth = 11- The effect of removing significant clus- 
ters is mostly apparent on intermediate scales. It is less 
pronounced than when cutting the pow er spectrum with 
respect to halo masses, as shown by iTakada fc Bridle! 
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(|2007D . because of the weak lensing selection function, 
that mixes masses when integrating the matter power 
spectrum on redshift. 

The power spectrum is a measure of the Gaussian 
field only, and one must rely on direct measures of non- 
Gaussianity to complement it, like the bispectrum or 
cluster counts. 

2.2.2. Bispectrum 

Three point statistics are the lowest statistics to cap- 
ture non-Gaussian features in a statistical field. The dark 
matter and the weak lensing bispectra have been t he ob- 
ject of numerous studies (e.g. iTakada fc JaiiJ2dM ). most 
commonly based on the perturbation theory on large 
scales and on the hyper-extend ed perturbation theory 
(HEPT) on smaller scales (e.g. iBernardeau et aLll2002l 
for a review). Here, we compare and combine the con- 
straints from the weak lensing power spectrum and bis- 
pectrum evaluated in the same halo model. This fills a 
gap in the literature, where such constraints come from 
the HEPT approximation only. 

The halo model allows us to separate the contribu- 
tion from the unresolved and the resolved parts of the 
field when computing the bispectrum, like we did for 
the power spectrum. However, discarding the unresolved 
part from the bispectrum removes too much useful sig- 
nal, and thus suppresses much power in the bispectrum. 
Therefore, we will always evaluate the bispectrum on the 
entire (resolved and unresolved) data field. 

2.2.3. Weak-lensing-selected halo counts 

Clusters of galaxies, tracing dark matter halos, have 
emerged from the primordial Gaussian field through 
non-linear gravitational clustering, and are those non- 
Gaussianities that one can detect with higher order 
statistics. Or one can simply estimate their abundance 
to take non-Gaussianities into account. In this paper, 
we concentrate on how to break the degeneracies from 
the power spectrum by combining it with non- Gaussian 
features. As a result, we will not focus on the constraints 
that halo counts alone can bring, but on how they im- 
prove the constraints from the power spectrum. 

Weak-lensing-selected clusters are detected on a mass 
map derived from the data used to measure the power 
spectrum and the bispectrum. Hence, they come at no 
extra cost, contrary to X-ray or optically-selected clus- 
ters. Defined as peaks with S/N greater than a cer- 
tain threshold vth, they are the resolved part introduced 
above. 

Once halos are selected, one can add information to 
their single S/N distribution, e.g. by introducing their 
redshift and mass, and thus measure their abundance as 
a function of those quantities. Hereafter, we investigate 
the constraints on cosmological parameters provided by 
four kinds of halo counts. With increasing information, 
we look at (1) counts as a function of halo's S/N only ; 
(2) counts as a function of halo's mass ; (3) counts as a 
function of halo's redshift ; and (4) counts as a function 
of halo's redshift and mass. 

In the remainder of this paper, we assume that the 
weak lensing selection function is perfectly known. We 
neglect the nuisance parameters commonly used to ac- 
count for the imperfect knowledge of the relation be- 



tween an observable and a cluster's mass (e.g. Im3[200l 
iLima fc Hul[200llMaiumdar fc Mohd[200l . 

2.3. Constraints on cosmological parameters 

Fisher matrices allow one to characterize how a set of 
observables x is able to constrain a set of parameters p 
around a fiducia l model. The associated Fisher matrix is 
defined by (e.g. iHu fc Tegma rk 1999) 



dpadpp 



(6) 



where L(x, p) is the associated likelihood. 

Given a fiducial model, the inverse of the Fisher ma- 
trix estimated in its neighborhood provides a lower 
limit to the parameters' covariance matrix. It thus quan- 
tifies the best statistical error tha t can be reached on the 
parameters, cr(pa) ^ \/ (F^^)acn where <j{pa) is the la 
error on parameter pa marginalized over other parame- 
ters pp. 

The Fisher matrices for the weak lensing power spec- 
trum, bispectrum and c l uster counts can b e found e.g. 
inlCoorav fc Hul (l20(nbh : iLima fc Hul (l200l:lHu fc Jain 
(120041) : ITakada fc JainI (|2004D and iTakada fc Bridle 
(l200l . 

We aim to determine how to best capture non- 
Gaussianities in order to break the degeneracies from 
the power spectrum alone. We thus need to combine the 
power spectrum with the bispectrum or the halo counts. 
To combine two observables Di and D2, we create the 
data vector D = {Di, D2}. The associated Fisher matrix 



is F, 



((l) + (2)) 



D 



The rigorous estimation of the combined observables' 
covariance matrix. 



C((l) + (2)) 



C(i) c(^)^C^) 
C(i)^(2) c(2) 



(7) 



requires the knowledge of the cross-covariance 

The weak lensing power spectrum is dominated by the 
most massive halos' contribution, those which are the 
most likely to be detected and taken into account in clus- 
ter counts. Hence, the power spectrum and halo counts 
are not independent observables : their cross-covariance 
is n on-zero, and mu s t be accounted for. Nonethe- 
less, ITakada fc Bridld ()2007l ) showed that neglecting it 
changes the errors by only a few percent. 

Here, we take a different approach, by separating the 
contribution of different halos to the power spectrum 
(Eq. ini). We combine the counts of clusters with S/N 
greater than the threshold vti, (the resolved part) with 
the power spectrum created by all halos with S/N smaller 
than the same threshold t'th (the unresolved part). If we 
neglect the clustering between halos, the power spectrum 
in the unresolved part and the resolved clusters are un- 
correlated : the covariance between cluster counts and 
the power spectrum vanishes, and we can simply add the 



Fisher matrices, f'^^^^^ 



afi ' 



where the su- 



perscript stands for the power spectrum evaluated 

on the unresolved part. 

As shown by ITakada fc JainI (|2004D . the cross- 
covariance between the power spectrum and the bispec- 
trum has no Gaussian feature, but arises from the five- 
point correlation function of the shear field. In such a 
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case, it is safe to approximate the Fisher matrix of the 
combination between the power spectrum and the bis- 
pectrum by the sum of the individual Fisher matrices, 

p(ps+bisp) ^ pps pbisp 

3. RESULTS 

In this section, we compute and compare the er- 
rors on cosmological parameters, that can be reached 
by measuring the power spectrum alone, and when 
combining it with either counts or the bispectrum. 
We assume a fiducial cosmology described by the 
parameters (ri™, ^a, ilb, <7s, h, n, wq, Wa) = 
(0.3, 0.7, 0.04, 0.9, 0.7, 1,-1, 0). The Universe's curvature 
is a free parameter, set by the combination of and 
n\. The evolution of the dark energy equation of state 
is parametrized by w{a) = wq + {1 — a)wa, a being the 
scale factor. We assume a Euclid-like survey, 20,000 deg^ 
wide, with ng = 40 galaxies arcmin"^ with median red- 
shift z„i = 1. We assume that the intrinsic shape r.m.s 
is CTint = 0.3. 

To count halos, we use 20 linear redshift bins, span- 
ning the entire interval accessible to our fiducial survey 
(0 < z < 5) - the highest bins being empty - and 15 
logarithmic mass bins (lO^Oh^^Mo < M < lO^^h-^M© 
- where the lower bound is small enough so that all clus- 
ters detectable by their weak lensing signal are accounted 
for) to compute halo counts as a function of redshift 
and/or mass, and 15 linear S/N bins for S/N peak counts 
(2 < S/N < 20). We estimate the power spectrum in 10 
redshift bins (0 < z < 5). For numerical reasons, we 
cannot use more than 3 redshift bins to compute the 
tomography bispectru m. This is not problematic, since 
iTakada fc JainI (|2004l ) showed that the tomography bis- 
pectrum S/N quickly converges to its maximum value, 
and almost reaches it for three redshift bins. 

Before re porting our resu lts, we would like to refer to 
the work of lVallisneril ([2008) in which he shows the limi- 
tations of the Fisher matrix formalism. In particular, one 
must pay particular attention to ill-conditioned Fisher 
matrices, the inversion of which is likely to be wrong. 
Here, we check the condition number of all of our Fisher 
matrices, and consider as good only those with a small 
enough condition number so that their inversion can be 
trusted. Wc find that the Fisher matrices for counts as 
a function of S/N, of mass and of redshift, are all ill- 
conditionned. Therefore, such halo counts cannot give 
reliable Fisher constraints by themselves. This is true 
for our particular set of parameters though, and could 
not be true on other parameter spaces. Looking for such 
spaces is beyond the scope of this paper, and we will give 
constraints from halo counts only when combined with 
the power spectrum, for which the condition number is 
low enough to be safe. 

The goal of this paper is to compare how cluster counts 
and the bispectrum capture non-Gaussianities and break 
degeneracies left over by the power spectrum. Conse- 
quently, we ignore combining cluster counts with the 
bispectrum, as well as combining the three observables 
(power spectrum, bispectrum and cluster counts) at a 
time. 

We evaluate Fisher matrices for power spectra esti- 
mated both on all the data and on the unresolved parts 
defined by different vth, for the corresponding halo counts 



(halos with v > t'th), and for the bispectrum, and the 
combinations introduced in section 12.31 Figure [3] shows 
the expected errors that they provide, marginalized over 
all eight parameters, as a function of the S/N thresh- 
old fth. The power spectrum Pk(< i^th) is estimated on 
the unresolved part (dashed red), then combined with 
clusters counts of clusters (made as a function of mass 
and redshift - solid green-, redshift only -solid blue-, mass 
only -solid cyan- and S /N only -solid purple) . The dotted 
black flat line shows the errors brought by the measure- 
ment of the power spectrum on the entire data set. Dis- 
carding the contribution of the most massive halos (i.e. 
the most non-Gaussian features) from the power spec- 
trum improves the errors on D,„i and erg i the errors being 
smallest for an optimal threshold I'th.opt ~ 6. Although 
the counter-intuitive increase in the performance that we 
observe when discarding the resolved part of the data is 
not so surprising, since the impact of non-Gaussiani ties 
on the power spectrum is lowered (jShaw et all 120091 ob- 
serve a similar trend for Sunyaev-Zel'dovich clusters), we 
should note here that it could be exacerbated if our anal- 
ysis is close to the Fisher matrices formalism's limits. 
Moreover, there may be an optimal weighting scheme 
that would both allow for the signal in the resolved part 
and minimize its non-Gaussian errors. Smaller errors 
on the cosmological parameters would then be expected, 
since more information (with minimal noise) would be 
taken into account than that we use here with the un- 
resolved part. Investigating these issues is beyond the 
scope of this paper, and we defer it to a later study. The 
dark energy equation of state parameters wq and Wa are 
less affected by the separation of the power spectrum into 
the resolved and the unresolved parts. 

Although cluster counts do not provide strong con- 
straints by themselves (see below. Table [T]) , combining 
them with the power spectrum improves the errors it 
provides, in particular when low S/N clusters are taken 
into account. When the clusters considered are signifi- 
cant enough, all four kinds of counts perform as well at 
capturing non-Gaussianities : the power spectrum pro- 
vides most of the Fisher information, and clusters help by 
breaking parameters degeneracies, independently of the 
way they are considered. Figure [3] shows that this is true 
for vth ^ 6, a saf e threshold to discard false detections 
(jPace et al.l l2007[ ). Therefore, that makes S/N peaks 
counts a direct and efficient prob e to combin e with t he 
power spectrum, as shown e.g. bv lPires et al.l (|2009bl ). 

The dash-dotted yellow line on Fig. |3] shows the er- 
rors given by the bispectrum, and the long-dashed orange 
line shows those given when combining it with the power 
spectrum. At low S/N thresholds, when combined with 
the power spectrum, the bispectrum and clusters counts 
give similar constraints on the four parameters ftm, fg, 
Wq and Wa- At higher thresholds the bispectrum gives 
better constraints than clusters counts on the dark en- 
ergy parameters, yet errors remain comparable for r2„j 
and as- 

Table [T] lists the marginalized errors on all of our eight 
parameters, as well as those of the dark energy equation 
of state Wp at the pivot point (section 13. 1|) . for the mea- 
surements and combinations considered, with vth = 6. 
Halo counts (excepted counts as a function of redshift 
and mass) are plagued by ill-conditioned Fisher matri- 
ces, and we do not consider them alone. Nevertheless, 
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Fig. 3. — Marginalized errors on cosmological parameters, as a function of the threshold between resolved and unresolved parts. Top-left: 
Qrn- Top-right : erg. Bottom-left ; uiq. Bottom-right : Wa- Different colors and line styles label different measurements : power spectrum 
on all the data (dotted black), power spectrum on the unresolved part (dashed red), bispectrum (dash-dotted yellow), and combinations of 
the power spectrum on the unresolved part with clusters counts as a function of redshift and mass (solid green), redshift only (solid blue), 
mass only (solid cyan), S/N only (solid purple) and with the bispectrum (long-dashed orange). Note that the y-scalc is not the same for 
top and bottom panels. 



TABLE 1 

Marginalized absolute errors on cosmological parameters. 











n(M,z) P, 


. " + n(M,z) 


P"''' + n(z) 


P,"'"' + n(M) P,"'"' + n(u) 


P^"" + B« 


Pk+B^ 


(0.3) 


0.0049 


0.0027 


0.0047 


0.070 


0.0024 


0.0024 


0.0026 


0.0026 


0.0021 


0.0023 


Ua (0.7) 


0.024 


0.012 


0.019 


0.254 


0.0099 


0.010 


0.011 


0.011 


0.0087 


0.0094 


Cfc (0.04) 


0.015 


0.015 


0.016 


0.656 


0.013 


0.014 


0.014 


0.014 


0.008 


0.008 


<T8 (0.9) 


0.0086 


0.0039 


0.0065 


0.112 


0.0032 


0.0034 


0.0036 


0.0036 


0.0029 


0.0033 


h (0.7) 


0.091 


0.086 


0.122 


3.59 


0.077 


0.079 


0.079 


0.085 


0.048 


0.051 


wo (-1) 


0.040 


0.041 


0.031 


0.283 


0.023 


0.024 


0.025 


0.025 


0.016 


0.016 


Wa (0) 


0.147 


0.133 


0.156 


2.179 


0.062 


0.063 


0.064 


0.064 


0.055 


0.063 


Wp 


0.022 


0.015 


0.013 


0.252 


0.014 


0.014 


0.015 


0.015 


0.0049 


0.0045 


nil) 


0.020 


0.021 


0.038 


1.00 


0.018 


0.018 


0.019 


0.020 


0.011 


0.012 



The S/N threshold between the resolved and unresolved parts is set to i-'th = 6. 

The central value for the parameter in our fiducial model is given into parenthesis in the first column. 

We ignore combining the power spectrum with the bispectrum and cluster counts at a time (see main text). 



since counts as a function of S/N, of mass alone and function of mass and redshift, their expected constraints 
redshift alone contain less information than counts as a are weaker than those listed in table [T] for counts as a 
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Fig. 4. — Dark energy FoM, as defined by Eq. as a function 
of the threshold between the resolved and the unresolved parts, for 
the same measurements as in Fig. |3] labeled with the same colors 
and line styles. Thick lines denote measurements made with the 
power spectrum estimated on the unresolved part, and thin lines 
take the entire information into account in the power spectrum 
estimation. 



function of mass and redshift. We also list in the ta- 
ble the errors provided when combining the bispectrum 
with the power spectrum estimated on all the data. Al- 
though this combination is not explicitly shown on figure 
[3l it can be drawn from the high S/N threshold-limit of 
the combination of the bispectrum with the power spec- 
trum estimated on the unresolved part (long-dashed or- 
ange line). The table emphasizes the conclusions given 
by figure [H The four combinations of clusters counts 
with the power spectrum give similar constraints, that 
are comparable to those provided by the combination of 
the bispectrum and the power spectrum. Furthermore, 
although the bispectrum and the power spectrum give 
comparable errors, combining them breaks degeneracies 
a nd lowers er r ors by a factor 2. 

iPires et "all ()2009b[ ) looked at how several statistics 
(clusters counts, bispectrum, skewness, kurtosis) break 
the flrn-o's degeneracy, and reported the high efficiency 
of cluster counts. They found that the skewness could 
also discriminate against models, but that the bispec- 
trum, when using only equilateral triangles, gave poor 
results. We expand their work by taking into account all 
triangle configuration when computing the bispectrum, 
and add it to clusters counts as an efficient way to break 
degeneracies (see section [3.2.2|) . In particular, Fig. [3] 
shows that this conclusion is true not only for the flm- 
CTg degeneracy, but also for the dark energy equation of 
state. 

3.1. Dark energy plane 

We now specialize the discussi on to the dark energ y 
figure of merit, as defined e.g. bv lAlbrecht et al.l (|2009l ). 



(T{Wp)<T{Wa) 



Fig. 5. — Constraints on the Qrn-crs plane, when all other param- 
eters are kept constant, for a 20,000 deg^ survey, and no tomog- 
raphy. Filled ellipses show constraints from the power spectrum 
(cyan), S/N peak counts (grey), cluster counts as a function of 
mass and redshift (green) and bispectrum when all triangle config- 
uration are taken into account (yellow); the constraints obtained 
from the bispectrum when taking into account equilateral trian- 
gles only get out of the frame. Open ellipses are constraints from 
combining the power spectrum with the bispectrum (all triangle 
configurations - blue), with peak counts (black) and cluster counts 
(red). 



where Wp is the dark energy equation of state at the pivot 
redshift, where the dark energy is best constrained. The 
errors on Wp are listed in Table [1] 

Figure m shows the dependence of this figure of merit as 
a function of the S/N threshold i/th between the resolved 
and the unresolved parts of the data. The color code 
and line styles are the same as those used for the previ- 
ous figure. Thick lines show the FoMs when the power 
spectrum is estimated on the unresolved part, and thin 
lines show the FoMs when it is estimated on all the data 
(we neglect the covariance between the power spectrum 
and cluster counts). 

Considering the power spectrum on the unresolved 
part only (dashed red) slightly improves the FoM. This 
holds true when combining clusters with > 5. The 
bispectrum gives a higher figure of merit than cluster 
counts, and depends very weakly on vth- Only at low 
S/N do the figures of merit compare. 

Clusters counts, when the selection function and the 
contamination are controlled well enough to allow one to 
consider a very low S/N threshold, and the bispectrum 
are comparable at extracting non-Gaussian aspects in 
a weak lensing survey and at constraining dark energy. 
On the other hand, if we use a safe threshold for clus- 
ters detection (z^th ~ 6), the bispectrum becomes better 
than clusters counts to break degeneracies and to mea- 
sure dark energy. 

3.2. Comparison with other works 

3.2.1. Peak counting 

In two recent papers , iDietrich fc Hartlad (|2009f ) and 
iKratochvil et al.l (|2009[ ) show with numerical simulations 
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how powerful weak lensing peaks counting is at con- 
straining and CTg, and w, respectively. These appear 
to contradict our results, in which we saw that if indeed, 
counting peaks efSciently captures non-Gaussianity and 
breaks power spectrum degeneracies, it does it only when 
combined with the power spectrum, but gives poor con- 
straints when used alone (Tabled]). A concern expressed 
by the aforementioned authors is that they varied only 
their parameters of interest (flm and cg or w), and defer 
marginalizing on other parameters for later works. In our 
analysis, we indeed marginalized on eight parameters. 
To show that this considerably lowers the constraints, 
we make the same kind of analysis, by varyin g $7^ and 
(78 only, as done by [Dietrich fc Hartlapl ()2009( ) . For sim- 
plicity, we do not separate the power spectrum into re- 
solved and unresolved parts, but consider the total power 
spectrum only, without tomography. Since our analytical 
model cannot fully reproduce their taking into account 
cosmologically significant false detections (peaks that do 
not correspond to a single, localized, virialized cluster, 
but that are instead due to mass alignments along the 
line of sight), we use counts as a function of S/N as a 
proxy of their measure. We use a detection threshold 
lyth = 3. Figure [S] shows how counting peaks (grey filled 
ellipse) compares to the power spectrum (cyan filled el- 
lipse) in that case. The black open ell ipse shows the com- 
binati on of the two. As shown by [Dietrich fc Hartlapl 
(|2009| ). peak counts indeed is an efficient cosmological 
tool when varying and erg only. Marginalizing on 
other parameters hampers its ability to perform, as we 
showed above, unless combined with the power spectrum. 

However, since our model d oes not catch filaments and 
signi ficant false detect i ons, a s [Dietrich fc Hartlapl (|2009[ ) 
and iKratochvil et al.l (|2009[ ) do, our conclusion about 
cluster counts alone may be over-pessimistic. Therefore, 
it will be valuable to adapt their simulations to higher 
dimensional parameter estimations. 

3.2.2. Bispectrum : the need for all triangle 
configurations 

Looking for the b est statistics to discriminate models, 
iPires etaLl (|2009bD conclude that peak counting is the 
best way of breaking the flm — cg degeneracy, with other 
third-order statistics (such the skewness) also being able 
to break degeneracies. In that work the bispectrum in 
particular was found to be weak at breaking degenera- 
cies. This apparent discrepancy with our findings here, 
where we find that combining the bispectrum with the 
power spectrum is comparable to combining peak count- 
ing with the power spectrum, comes from the fact that 
the bispectrum presented here uses all tria ngle configu- 
rations, while the bispectrum in lPires et al.l (|2009bj used 
only equilateral triangles. 

To emphasize the importance of considering all tri- 
angle configurations when computing the bispectrum, 
we compute the constraints expected on and ag, 
while k eeping all other parameters fixed, as in lPires et al.l 
(|2009bD . from the power spectrum, the bispectrum with 
all triangles information or only equilateral triangles, 
and cluster counts. We use one redshift bi n. Our an- 
alytical model does not allow us to compute iPires et al.l 
(|2009h) 's Wavelet Peak Counts (counts of clusters in dif- 
ferent wavelet scales), so we use the counts as a function 
of mass and redshift as a proxy. 



Figure |5| shows how constraints from these measure- 
ments compare. The cyan filled ellipse shows constraints 
from the power spectrum, the green one, constraints from 
cluster counts, an d the red open ellipse their combina- 
tion. As shown by IPires eflil. (2009b), cluster counts are 
very efficient at grabbing non-Gaussianities and breaking 
the power spectrum degeneracies in that particular case. 
The yellow filled ellipse shows constraints from the bis- 
pectrum, when the information from all triangle config- 
urations is taken into account, and the blue open ellipse 
shows constraints from combining this bispectrum with 
the power spectrum. As mentioned above, the bispec- 
trum (with information from all triangle configurations) 
provides tight constraints, and is as efficient as cluster 
counts at capturing non-Gaussian features. However, we 
find that when taking into account only equilateral tri- 
angles when computing constraints from the bispectrum, 
the errors on the parameters reach 100%, making such 
a measurement extremely low efficiency. The frame of 
Figure \5\ is much too small to sho w the associated el - 
lipse. Consequently, we concur with IPires et~all (|2009bD 
when claiming that the bispectrum measured with equi- 
lateral triangles only does not allow one to break the 
power spectrum degeneracy between rim and erg . 

Coming back to our 8-parameters analysis, we notice 
the same behavior as that mentioned in the above para- 
graph. To capture non-Gaussian features and efficiently 
break degeneracies left over by the power spectrum, one 
has to take into account all triangle configurations when 
measuring the weak lensing bispectrum. Equilateral tri- 
angles only do not provide any cosmological constraints. 

4. CONCLUSION 

We investigated how to best capture non-Gaussianities 
in weak lensing surveys so as to break degeneracies left 
over by the power spectrum in weak lensing analysis, by 
means of combining it with galaxy clusters counts and 
with the weak lensing bispectrum. To this end, we com- 
puted, in the halo model framework, the power spectrum, 
the bispectrum and four kinds of halo counts (as a func- 
tion of S/N, mass, redshift, and mass and redshift), in a 
Euclid-like survey. We have computed the errors on cos- 
mological parameters brought by the combination of the 
power spectrum with the bispectrum, and of the power 
spectrum with clusters counts. We have emphasized our 
discussion on the varying dark energy equation of state, 
the mass density parameter 0^ and the amplitude of 
density fluctuations erg, marginalized on a 8-parameter 
cosmological model. We did not take any prior from 
outside weak lensing (e.g. CMB priors), so as to show 
how a weak lensing survey alone can allow us to capture 
non-Gaussianities, and what constraints we can expect 
from it on cosmological parameters. 

When computing the power spectrum, we have sepa- 
rated the contribution of the resolved and the unresolved 
parts of the data, delimited by a S/N threshold i^th. Do- 
ing so, we noted an improvement in errors on cosmolog- 
ical parameters. 

Although peak counts alone are powerful at constrain- 
ing f2„j and erg wh en other parameters are kept con- 
stant (as shown by [Dietrich fc H artlap 2009^). we have 
seen that marginalizing on a 8-parameter space make 
them provide only weak constraints, unless combined 
with the power spectrum. We have indeed shown that 
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when combined with the power spectrum, cluster counts 
provide percent level marginalized errors on cosmologi- 
cal parameters. Furthermore, provided that the thresh- 
old is high enough (i^th ~ 6), all four kinds of cluster 
counts we consider give similar constraints. Hence, com- 
bining the power spectrum with S/N peak counts gives 
constraints nearly as tight as combining the power spec- 
trum with clusters counted as a function of redshift and 
mass. Galaxy clusters provide the same non-Gaussian in- 
formation, independently of the way they are taken into 
account, as long as they are combined with the power 
spectrum : different cluster counts break degeneracies in 
a similar way. This will allow us to put stringent con- 
straints on cosmological parameters from weak lensing 
alone, by combining S/N peaks with the power spectrum, 
without any extra cost and requirement of follow-up for 
clusters' redshift and mass measurement. 

We have shown that the weak lensing bispectrum ef- 
ficiently captures non-Gaussianities, and is marginally 
more efficient than cluster counts at breaking degen- 
eracies left over by the power spectrum. Furthermore, 
the bispectrum alone provides constraints comparable to 
those giveii by the power spectrum. We have expanded 
iPires et aD (|2009bl )'s conclusions by adding the bispec- 
trum to clusters counts as efficient tools to capture non- 
Gaussianities. We also noted that equilateral triangles 
alone do not bring significant information. To provide 
percent level constraints, all triangle configurations must 
be accounted for when measuring the weak lensing bis- 
pectrum. 

All our results are based on a Euclid-like survey, and 
may not be directly applicable to current surveys such 
as the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) or the Cosmic Evolution Survey (COSMOS). 
In particular, on current surveys, the cosmic variance 
could play an important role in lowering the ability 
of the bispectrum to be more efficient than S/N peak 
counts. How those two measures compare to capture 
non-Gaussian features could then depend on the survey's 
characteristics. This question is beyond the scope of this 
paper, and will be addressed in a subsequent work. 

We have assumed that all systematics were well ac- 
counted for so as not to bias our results. Although that 



can seem over-optimistic, the amount of efforts currently 
underway to correct for various systematics gives us good 
reasons to think that our assumption is likely to be met 
when a Euclid-like survey is undertaken. Moreover, al- 
lowing for systematics requires a more sophisticate d anal- 
ysis than the Fisher matrix (jAmara fc Refregied 12008). 
We defer it to a future paper, looking at how systematics 
enter in the correlations between different statistics. 

In light of our results, it appears that weak lensing sur- 
veys alone will be able to reach the percent accuracy on 
the dark energy equation of state parameters. Moreover, 
it provides us with ways to cross-check the parameters' 
measurements. For instance, after measuring the power 
spectrum, counting S /N peaks provides an easy and fast 
way to optimally capture non-Gaussianities and tighten 
constraints. Then, the bispectrum and its combination 
with the power spectrum should give consistent parame- 
ters' estimation and similar constraints. However, mea- 
suring the bispectrum on half the sky will constitute a 
real challenge. Along with already existing robust ari d 
fast algorithms parvis et al.ll2004t IZhang fc Penll2005D . 
new ones are being developed (e.g. Berge et al in prep 
and Semboloni et al in prep) to measure either the 3PCF 
in real sp ace or the b ispectrum in Fourier space (as pro- 
posed by Fires et al. 2009a). In the meantime, more and 
more precise numerical simulations will be a needed tool 
to test and tune those algorithms, and current surveys 
such as the CFHTLS and COSMOS are ideal benchmark 
to test codes in the real world. 
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APPENDIX 
A - HALO MODEL 

Presented as an alternative to fitting functions bas ed on numerical simula tions, the halo model for large-scale 
structures describes structures as spherical haloes - see ICoorav fc ShethI (|2002r ) for a review. In this framework, the 
large-scale structures statistics can be described by correlations within a same halo (the profile of which affects them) 
and between different halos (the clustering of halos is thus naturally taken into account). Although being simplistic, 
the model agrees fairly well with common numerical simulations. Here, we summarize its key ingredients : the mass 
function, the profile of halos and the halo bias. We also give the expressions for the 3D dark matter power spectrum 
and bispectrum. 

Mass function 

The comoving number of virialized halos, whose mass ranges between M and M + dM and redshift ranges between 
z and z + dz, is given by (jPress fc Schechteilll974[ ) : 

niM,z)dM^^'^^mdM (Al) 

where po is the Universe's current background density and J'(M, z) — 6c{z)/a{M). 5c{z) is the non-linear overdensity 
of a halo collapsing at redshift z. In a dar k energy universe, Sdz) depends weakly on cosmology, especially on the 
dark energy's equation of state w. We use [Weinberg fc Kamionkowskil ()2003l )'s fitting function to compute Sc{z;w). 
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This fitting function is valid for constant w only, therefore we compute it using an effective equation of state WcS (a) = 
wq + {1 — a)wa- The r.m.s of the mass fluctuations in a sphere containing a mass M at redshift z is defined by 

where P''"(fc) is the three-dimension linear matter power spectrum, R — (3Af/47rp)^/'^ is the radius of the considered 

sphere, and p is the mean background density. The window function VFfx) = (3/x ) [sin(a:) — xcosjx)]. 

A couple of mass functions have be en described in t he literature (IPress fc SchechteiillQTl iSheth fc Tor meni 119991 : 
iJenkins et a"ni2001l : iTinker et all 120081 ) . We use that of iJenkins eFall (|2001[ ). f(a) = exp(-| Ina'^ + bjl"^), where 
aj = 0.315, bj = 0.61 and Cj = 3.8. 

Profile 

The top-hat collapse model for structure formation describes the matter infall on the gravitational wells until 
virialization. At this stage, non-linear physics t akes place, tha t furth er modifies the virialized halo profile, and that 
must be assessed through numerical simulation. iNavarro et al.l ()1996D showed that the mass profile of such virialized 
halos can be described by a "universal" profile, 

'"('^1^) = , / I (A3) 

where and Ps correspond to a characteristic radius and density, respectively, and (a, /3)=(1,2) for the usual NFW 
profile, that we use in this paper. 

The mass of a NFW halo is then given by M = dr47rr^p(r|M), where rvir is the halo's virial radius. In 
theory, a NFW is infinite, but we assume here that it is truncated at rvir- The virial and characteristic radiu s are 
linked by ry ir = ct s , where c is the halo's concentration, that can be parametrized, following iBuUock et al.l (|2001h and 
iDolag et all (|2004[ ). by 

c(M, z) = ^( ^^^'^ . (A4) 

^ ' ' l + z \M,{z = 0)J ^ ' 

In this equation, M h.(z = 0) is the non-linear mass scale, defined as i^(M, ,z = 0) = 1. In contrast to earlie r work s 
jBullock et al."2001), wc a dd an e xtra redshift-dependence in the parameter /3. Based on IZhao et al.l (|2003l . [2"008[) : 
[DufFy ct al. (2008) and Ga o et ahl ([^08), it allows us to avoi d the catastrophic dr op in concentration observed on 
simulations for high redshift, massive halos, as parame trized by [Bullock et al.l (120011) with a c onstant /?. We find that 
Co = 8 and P(z) = 0.13- 0.05z, besides agreeing with lDuffv et al.l (|2008f ) : iGao et al.l (|2008l ): IZhao et all (I2008D. and 



with iBullock et all ()200lD at low redshift, gives us the best agreement between the halo model and the ISmith et al.l 
(j2003t )'s power spectrum and bispectrum. One should note that the concentration has a significant impact on the 
1-halo term of the halo model power spectrum and bispectrum, as shown by .Huffcnberger fc Seljak (,2003) . and thus 
a thorough description is needed. 

Bias 

The bias describes how ha los cluster with respect to the matter distribution. It is well described at first order by 
fitting functions of the form ()Mo fc Whitelll996HSheth et al.ll200H) . 

b{M,z) = 1- —— + ——— (A5 
6c[z) 1 -I- {qiy)P 

with q = 0.707 and p — 0.3. We neglect the higher order terms of the bias. 

Matter power spectrum and bispectrum 

The halo model describes the correlation functions of the density field, in real-space, as the sum of the correlation 
between points belonging to a same halo and of the correlation between points in different halos. The same contributions 
appear in Fourier space, when defining the power spectrum a nd bispectruni, which are thus the sum of a 1-halo term 
and a 2-halo term (and 3-halo term for the bispectrum) (e.g. ICoorav fc Hull2001aD . Ps{k) = ''(fc) + Pg'^ik). 

The 1-halo component is given by 



Ps'^ik) = J dmn{m) (^j 



2 

2 



{k\m)\' (A6) 



and the 2-halo term is given by 



Ps^i^) = j dmin(mi) \ u{k\mi) j dm2n(m2) \ u{k\m2)Phh{k\mim2). (A7) 

The function 

, , > f^^" nSinfcrp(r|m) , 

u(km)= drAwr^— ' ' A8 

In kr m 
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is the Fourier transform of the dark matter density profile. The term Phh{k\'in'ifn2) represents the power spectrum 
of halos of mass mi and m2 and is approximated by Pfiii{k\mim2) ~ n^^-^&i(TOi)P]™(fc), where bi{mi) is halo i's bias 
and P]'"(fc) is the matter linear power spectrum. We use the lEisenstein fc Hul ()1998[ )'s transfert function to evaluate 

The 3D matter bispectrum can be decomposed into terms coming from one, two and three halos : 



(A9) 

(AlO) 
(All) 

(A12) 

where CP. denotes cir cular permutations. The J function is given by the second-order perturbation theory 
(jBernardeau et al.ll2002f ) : 



with 



and 



Bl''{ki,k2,h)=I^{ki,k2,k3), 

BfikiMM = ll{k^,k2)ll{k3)P,{k:i)+C.P., 

Bf{ki,k2,k3) = [2J{ki,k2,k3)ll{k3)+lf{k3)] X ll{ki)ll{k2)Ps{ki)Ps{k2)+C.P, 



J(fcl,fc2,fc3) = l-^f^™'/'' + 



7,2 „ 1.2 „ ^,2^ 2 
'^3 '^1 "'2 

2fcifc2 



kf 



1.2 _ )l2 _ 1,2 ' 7 
3 1 2 



(A13) 



The function is defined through the Fourier transform of the halo profile u{k\ni) (Eq. IA8I 



I^{ki, . . . ,k^; z) = I dm ^ — j n{m, z)bp{'m) x u{ki,m) . . .u{kf^,'m) 



(AM) 



with 6o = 1 



B - SIGNAL-TO-NOISE RATIO FOR A NEW HALO 

Here, we specialize the discussion of section 12.11 to an NFW halo. To evaluate equation ([1]) in this case, we first 
need to derive the convergence of such a halo. The mass density, projected along the line-of-sight, of a NFW halo at 
redshift z, with concentration c, is given by : 



dzp{x,z) = 2psrsg{x) 



(Bl) 



where x = r/rg- The func tion g depends only on the distance to the halo's center, and is given by (as long as c ^ 1) 
(|Wright fc Brainerdl[2000l) 



'j{x) = < 



(l-x^)(l+c) ' (l-2;2)-V 
3(l+c) \, ^ 1+c 



rarccosh 



x{l+c) ' 
(X< 1) 



(X=l) 



I 



(1 < a; ^ c) 
{x > c). 



The halo's convergence is then, if we assume that all sources are at the same redshift Zg 



Zs) = ^ = 2- g{x), 



where the critical density is defined as 



Scrit — 



crit ^crit 

AttG DdDds ■ 



(B2) 



(B3) 



(B4) 



In equation (jB4[) . c is the speed of light, G the gravitation constant, and Dg, and D^s are the angular-diameter 
distances of the source, of the lens, and between the lens and the source, respectively. 

To take the redsh ift distributio n Pzizs) of source galaxies, wc follow iSeitz fc SchneidCTl (flOOl . 

iBartelmann &: Schneideil ()2001|) and iWeinberg fc Kamionkowskil (|2002D and introduce the function 



Z{zs;zd) 



Hm, 



3 ^crit [Zd] Zs) 



-'crit,cxD 



Scrit(-Zdj Zs) 



l{Zd]Zs)' 



(B5) 



where Zd is the lens' redshift, and Scrit,oo is the critical density for a source at infinite redshift. This function 
allows us to project sources with a certain redshift distribution on a single redshift Zs verifying Z(zg) = (Z), with 
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(Z) — J (lzsPz{zs)Z{zs] Zd). One can then link the halo's convergence with a source at Zg to that with infinite redshift 
source, k{x, Zg) = k(x)Z(zs', Zd) (jSeitz fc Schneideilll997t) . The convergence given by a source population distributed 
in redshift along Pz{zs) is then given by k{x,Pz{zs) = (Z) k{x). 
The signal-to- noise ratio (Eq. [T|), for a NFW halo is then 

z^NFW = 2^/2^ (Z) ^ jj^^ y/G{c) (B6) 

where G{c) = Jq dx xg{x)^ is a function of the halo's concentration only, and is well fitted, for 1 ^ c ^ 200 (the 
maximum deviation being of order 0.5%), by 

0.131 0.375 , 7 , , , 

G c) w -I- 0.388 - 5 X 10^-^c - 2.8 x lO^V. (B7 

c 
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